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We consider a smooth forward facing step defined by the Gauss error function of height 
4-30% and four times the width of the local boundary layer thickness (^ 99 . The boundary 
layer flow over a smooth forward-facing stepped plate is studied with particular em¬ 
phasis on stabilisation and destabilisation of the Tollmien-Schlichting (TS) waves and 
subsequently on transition. The interaction between TS waves at a range of frequencies 
and a base flow over a single/two forward facing smooth steps is conducted by linear 
analysis. The results indicate that for a high frequency TS wave, the amplitude of the 
TS wave is attenuated in the unstable regime of the neutral stability curve correspond¬ 
ing to a flat plate boundary layer. Furthermore, it is observed that two smooth forward 
facing steps lead to a more acute reduction of the amplitude of the TS wave. When the 
height of a step is increased to more than 20 % of the local boundary layer thickness for 
a fixed width parameter, the TS wave is amplified and thereby a destabilisation effect is 
introduced. Therefore, stabilisation or destabilisation effect of a smooth step is typically 
dependent on its shape parameters. To validate the results of the linear stability analysis, 
where a high-frequency TS wave is damped by the forward facing smooth steps direct nu¬ 
merical simulation (DNS) is performed. The results of the DNS correlate favorably with 
the linear analysis and show that for the investigated high frequency TS wave, the K- 
type transition process is altered whereas the onset of the H-type transition is postponed. 
The results of the DNS suggest that for a high-frequency perturbation = 150 and in 
the absence of other external perturbations, two forward facing steps of height 5% and 
12% of the boundary layer thickness delayed H-type transition scenario and completely 
suppresses it for the K-type transition. 
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1. Introduction 

1.1. Motivation behind the study of steps in boundary layers 

In environments with low levels of disturbances, transition to turbulence is initiated 
by the exponential amplification of the Tollmien-Schlichting (TS) waves followed by 
the growth of secondary instabilities. Breakdown to turbulence generally occurs when 
the amplitude of the primary instability typically reaches 10 % of the free-stream veloc¬ 


ity magnitude (Herbert 1988 Cossu & Brandt 2002 ). The classical process of laminar- 


turbulent transition is subdivided into three stages: receptivity, linear eigenmode growth 
and non-linear breakdown to turbulence. A long-standing goal of laminar flow control 
(LFC) is the development of drag-reduction mechanisms by delaying the onset of tran¬ 
sition. The process of laminar to turbulent transition has been shown to be influenced 
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by many factors such as surface roughness elements, slits, surface waviness and steps. 
These surface imperfections can significantly influence the the laminar-turbulent transi¬ 
tion by influencing the growth of TS waves in accordance with linear stability theory and 
then non-linear breakdown along with three dimensional effects ( KachanQ^|1994 ). Since 
the existence of TS waves was confirmed by Schubauer & Skramstad (1948), numerous 


studies aiming to stabilise or destabilise the TS modes have been carried out in order 
to explore and explain different paths to transition. If the growth of the TS waves is 
reduced or completely suppressed and, providing no other instability mechanism comes 


into play, it has been suggested transition could be postponed or even eliminated (Davies 


& Carpenter 1996). Despite roughness elements being traditionally seen as an impedi¬ 


ment to the stability of the flat plate boundary layer recent research has shown this 


might not always be the case. For example, Shahinfar et ai (2012) showed that classical 


vortex generators, known for their efficiency in delaying, or even inhibiting, boundary 
layer separation can be equally effective in delaying transition. An array of miniature 
vortex generator (MVGs) was shown experimentally to strongly damp TS waves at a 
frequency T = x 10^ = 102 and delay the onset of transition. Similar results 

were obtained for T =135 and 178. Downs & Fransson (2014) found that TS wave 


T G {100,110,120,130} amplitudes over spanwise periodic surface patterns can be re¬ 
duced and demonstrated substantial delays in the onset of transition when TS waves are 
forced with large amplitudes. 

Over the past two decades, most investigations on topics of laminar-turbulent transi¬ 
tion in a boundary layer have focused on two kinds of problems: the receptivity mecha¬ 
nism ( |Wu|2001a|&l [Saric et a/.|2002||Ruban et a/.]20i^ a nd stabilisation/destabilisation 


of TS waves ( ^ossu fc Brandt [2 002 2004 Fransson et a/.|[^05[ 2006[ [Garzon fc Roberts] 

2013). Receptivity is the initial stage of the natural transition process, first highlighted by 


Morkovin (1969a), where environmental disturbances, such as acoustic waves or vortic- 


ity, are transformed into small scale perturbations within the boundary layer (Morkovin 


19696). The aim of these studies is to assess the initial condition of the disturbance am¬ 


plitude, frequency, and phase within the boundary (Morkovin 19696 Saric et al. 2002) 


So far, the receptivity mechanism of isolated small height roughness is well understood 
from theoretical, numerical and experimental points of view( Gaster|1965 Murdock|1980 


Goldstein|[T^ |Kerschen|[T9^ |W?I||20016| [Saric et a/.||2002D . The re- 


ceptivity mechanism shows that the deviation on the length scale of eigenmodes from 
a smooth surface can excite TS waves by interacting with free-stream disturbances or 


acoustic noise. From a theoretical point of view, Ruban (1984); Goldstein (1985); Duck 


et al. (1996) studied the interactions of free-stream disturbances with an isolated steady 


hump within the viscous sublayer of a triple-deck region. However, theoretical studies 
of the interaction between the TS wave and a distorted base flow have received less 
attention. More recently, the important theoretical work on the interaction of isolated 
roughness with either acoustic or vortical freestream disturbances was investigated by 
Wu & Hogg ( 2QQ6|) with in the framework of the triple-deck theory. For distributed rough¬ 


ness 


Corke et al. (1986) further inferred that the faster growth of TS waves on the rough 


wall was not attributable to the destabilization effect of roughness, such as an inflectional 
instability, but claimed that the growth was due to the continual excitation of TS waves 
on rough wall by free-stream turbulence. As we know, from a theoretical point of view, 
the understanding of the effect of isolated localised roughness on receptivity mechanism 
is better than that on growth of the TS wave ( Wu||2QQra Wu fc Hogg||2QQ6 ). 

In this paper, we investigate the effect of a smooth forward facing step on the growth 
properties of TS wave excited by forcing the boundary layer at different unstable non- 
dimensional frequencies. The amplitude of the forcing was chosen such that the velocity 
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profiles of the resulting TS waves are well resolved but weak enough that no secondary 
instabilities are introduced. For the domain of the fiat plate considered in the test cases 
presented here the unstable frequencies span IF G [27, 250] for displacement thickness 
Reynolds number Res* G [320,1500]. We denote frequencies as high for IF G [100,250]. 
These high-frequencies are of particular interest because they can lead to transition 


towards the leading edge of the fiat plate (Downs & Fransson 2014). Note, the differ¬ 
ence between the smooth step considered in this paper and a traditional sharp step 
is that the geometry for a smooth step is described by no less than two parameters 
(height and width) whereas for a traditional sharp step only one parameter (height) is 
required. |Nenni fc Gluyas (1966) first explicitly gave a critical height for sharp forward¬ 
facing steps corresponding to a Reynolds number, defined according to the step height. 


of Re/z^criticai = 1800. As discusscd by Edelmann & Rist (2015), the research on the influ¬ 
ence of steps on the stability of the boundary layer can be divided into two approaches: 
one focuses on finding a critical step height Recritical ( Nenni fc Gluyas]|1966 ) whereas 
the other is bas ed on the idea that the effect of a protuberance can be incorporated in 
the e^ method (Perraud et a/. [2004 Wang fc Gaster|20'Q5 Crouch et a/. [2006 Edelmann 


fc Rist|[2Q13 2015). It is worth mentioning that Wu fc Hogg (2006) showed that as the 

TS wave propagates through and is scattered by the mean-flow distortion induced by the 
roughness, it acquires a different amplitude down stream. They introduced the concept 
of a transmission coefficient. A further numerical study by Xu et al. (2016) confirms 


the localised isolated roughness has a local stabilising effect but overall a destabilising 
effect. Additionally an alternate expression of the transmission coefficient is introduced 
which can be incorporated into the e^ method. Recently, based on the second approach. 


Edelmann fc Rist (2015) found that generally, for transonic flows, sharp forward-facing 


steps led to an enhanced amplification of disturbances. They also found sub- and super¬ 
sonic results showed significant differences in the generation mechanism of the separation 
bubbles. A different phenomenon was found for incompressible flows when investigated 
numerically, by using an immersed boundary technique ( Worner et a/.|2003 ). The authors 
observed that for a nondimensional frequency R = 49.34, the amplitude of the TS wave 
is reduced throughout the domain considered by the forward-facing step. They attribute 
this stabilising effect to the thinner boundary layer evolving on the step in comparison 
to the boundary layer without a step. They also claimed that when a small separation 
zone appears in front of the step it has no influence on the TS wave. 

The presence of separation bubbles gives rise to a destabilising effect on a boundary 
layer. Generally, the separated shear layer will undergo rapid transition to turbulence 
and, even at rather small Reynolds numbers, separation provokes an increase in velocity 
perturbations and laminar flow breakdown, taking place in the separation region or close 
to it. The first observation of laminar separation bubbles were done by 
and the structure of a time-averaged bubble was given by 


Horton 


Jones (1938) 


(1968) and the in¬ 


terested reader can And a detailed review of the experimental work on the subject in 


Young fc Horton (1966). Hammond fc Redekopp (1998) found that a separation bubble 


could become absolutely unstable for peak reversed flow velocity in excess of 30% of 


the free-stream velocity. Theofllis (2000) report the shape of the globally unstable mode 
in recirculation bubble. A separation bubble can have important impact on the global 
stability of a boundary layer. In the following, we shall only focus on amplification of 
the TS wave by a separation bubble. Numerically, with laminar separation bubbles, |Rist| 
(1993) suggested a three-dimensional oblique mode breakdown rather than a secondary 


instability of finite-amplitude two-dimensional waves. Xu et al. (2016) investigated the 
behavior of TS waves undergoing small-scale localised distortions and found even a small 
separation bubble can amplify a TS wave. When a sharp forward-facing step of sufficient 
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height is present in a boundary layer, a separation bubble can easily be generated. In 
particular, the effective transformation or scattering of the freestream disturbances to 
the TS waves occurs preferably alongside sudden changes of the mean ffow (e.g. over the 
leading edge, separation region or suction slits). The external acoustics, surface vibra¬ 
tions and vortical disturbances in the form of localised ffow modulations or freestream 
turbulence are those which most frequently contribute to the boundary layer receptiv¬ 
ity as reviewed by Nishioka & Morkovin (1986); Kozlov & Ryzhov (1990); Saric (1990); 


Bippes (1999). Additionally, a smooth step is less receptive than a sharp step (Kachanov 
et a/.||1979 ). Another benefit of using a smooth step is to circumvent biglobal instability 


^ Hammond fc Redekopp|p^M ). 


1.2. Problem definition 

The smooth step considered in this paper, shown in figure is located in the unstable 
regime of the neutral stability curve close to the leading edge. This work mostly considers 
medium and high frequencies E G {100,140,150,160} TS waves but also the case of a 


low frequency forcing E = 49.34 is investigated for closer comparison with Worner et al 


(2003). The low frequency forcing is particularly interesting because it offers, in the 
context of the neutral stability curve of the zero pressure gradient flat plate, a much 
larger unstable regime compatible with the so called critical amplification factor N = S 
( Edelmann fc Ri^|2015 ). Recently, Downs & Fransson ( |2014| ) studied TS wave growth 
over spanwise-periodic surface patterns excited at E G {100,110,120,130}. They report 
that TS waves excited by high frequencies and large amplitudes < 0.48%l/oo): 

producing well resolved profiles without triggering secondary instabilities, can be reduced 
by spanwise-periodic surface patterns compared to the fiat plate case. Therefore, over 
a smooth step, understanding growth properties of the TS waves with high frequencies 
is pertinent. The linear analysis shows that in the presence of a single smooth step 
the TS wave can be attenuated below a critical height and smoothness but amplified 
above this critical height. Furthermore, the linear stability investigation of two isolated 
smooth steps instead of a single step with the same geometrical configuration revealed 
further reduction of the amplitude of the TS wave can be obtained compared to a single 
smooth step when the same step parameters bar location are used. Again, past a critical 
height, two forward facing smooth steps can have a destabilising effect on the TS mode. 
Furthermore, a smooth step can tolerate a higher height scale compared with a sharp step 
and does not introduce separation bubbles. In order to validate the stabilising effect of 
the smooth forward facing steps seen in the linear stability analysis, full non-linear direct 
numerical simulations of both the K- and H-type transition scenarios are conducted. 
These two transition scenarios are particularly compelling because they exhibit a long 
region of linear growth particularly suitable for the investigation of the effect of the 
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smooth forward facing step in the boundary layer on the TS wave ( Sayadi et a/.||2Ql'3 ). 
The results from the direct numerical simulation confirm the findings from the linear 
analysis. For both K- and H-type transition scenarios the forward facing smooth steps 
configuration has a stabilising effect even avoiding transition for the K-type scenario. 

The paper is organised as follows. In §2, we give fundamental definitions and describe 
the numerical tools employed in §3.1. In §3.2 and §3.3 and 3.4 we present results of linear 
analysis for high, medium and low frequency perturbations. §3.4, the results from the DNS 
are presented and discussed. A further discussion is then given in §4 and subsequently, 
the work is concluded. 


2. Mathematical formulations 

2.1. Fully nonlinear and linearised Navier-Stokes equations 

The non-dimensional momentum and continuity equations governing unsteady viscous 
flow with constant density are given as follows 

... “ " ( 2 . 1 ) 


dtUi - Re ^djUi + UjdjUi + diP = 


djUj = 


0 , 

0 , 


where Ui is one component of velocity field along the ith direction, dt denotes the 
derivative with respect to time, dj is the jth direction spatial derivative , Re is the 
Reynolds number defined by LUoq/u where L is the distance from the leading edge, v 
is the kinematic viscosity and p is the pressure. For a two-dimensional (2D) problem, 
i = 1,2 and ( 1 ^ 1 , 142 ) = {u^v) and for a three-dimensional (3D) problem, i = 1,2,3 and 

(U1,U2, W 3 ) = {U,V,W). 

Considering a steady state Ui of (13 about which a small perturbation such that 
Ui = UiFUi and dropping the second order terms in 12^, (2.1) can be linearised as follows 


dtUi — Re ^d‘iUi-\- UjdjUi + UjdjUi + dip = 


^3^3 ~ 


( 2 . 2 ) 


With suitable boundary conditions, in a linear regime, the system (2.2) can be used to 


exactly simulate evolution of a small perturbation Ui in a boundary layer. 

In the flat-plate simulations undertaken, with the assumptions of relatively large Re 
and no pressure gradient the base flow can be approximated by the well-known Blasius 
equation 


fiv) + ^f{v)f”{v) = 0 , 

subject to the following boundary conditions 

fiv) = fiv) = 0 at r? = 0, /' = 1 at r? 


(2.3) 


(2.4) 


where the prime denotes the derivative with respect to the similarity variable rj. Specifi¬ 
cally in the above, the dimensionless variables are defined by 


/ = T /\/nUooX and r] = y/\/Uoo/{i^x), 


(2.5) 


where T is the stream function. The streamwise and vertical velocity profiles of the 
Blasius boundary layer can be calculated by 

ub = Uoof'irt) and vb = ^]/~ /(^))- (2-6) 

Under the assumption of streamwise parallel flow in two dimensions, the perturbation 
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assumes the normal form 


{u, v,p) = {u, v,p) exp{i{ax — ujt)) + c.c., 


(2.7) 


where a and uj denote wave-number and frequency of a perturbation, respectively. The 
mode {u, v,p) in (2.7) generally can be obtained by solving the well-known Orr-Sommerfeld 
(0-S) equations, the solution of which for eigenvalues and eigenfunctions has been well 
studied ( |Stuart||1963 Schlichting fc Gersten|196^ Drazin & Reid'1981). When a surface 
imperfection occurs, the same notations u and v are used to denote the TS mode. 

Generally, for an unstable frequency uj G assuming that the TS mode is dependent 
on both X and y, the TS wave envelope is defined by the absolute maximum amplitude 
of the TS wave as follows 


A{x) = max{|{i(x, 7 ^,t)| : G [0,oc),Vt G 


( 2 . 8 ) 


2 . 2 . Definitions correlated with a surface imperfection 
In order to rescale the step, we introduce a reference boundary layer thickness Sgg = 
A.dlXcRef^^ and displacement thickness scales = 1.720SxcRefy‘^, defined accord¬ 
ing to a fiat plate boundary layer, where Xc is the distance from the leading edge to the 
centre position of a surface imperfection and Rcx^ = t/oo^c/^- We also let Res* = Uoo^* 1^ 
be the displacement Reynolds number. Now, we consider a forward-step-like surface im¬ 
perfection, which is defined by 

/.(.Y,i) = ^.(l + erf(^)) (2.9) 

where d and h{> 0 ) are the streamwise width scale and the normal direction length scale 
defined by the corresponding physical scales d and h as 

h = h/6gg, d = d/6gg, (2.10) 

and X is a streamwise local coordinate defined as follows 


X = {X- Xc)/6gg. 


( 2 . 11 ) 


For multiple smooth steps, the wall profile is formally defined by 

n 

Y,fs{X-Xi,h), ( 2 . 12 ) 

i=0 


where Xi {Xq = 0) denotes the relative centre position of each individual step with 
respect to the first step and n + 1 is the number of steps. 

In order to characterise or quantify geometrical steepness of a continuous function 
f{x), assuming at least f{x) G we introduce the following quantity 


j{x) = 


maxima;/(a;) I 


ymaV|a)(70r)p~+e^ 


-2 ^ 


(2.13) 


where 5 is a smooth parameter. It is clear that 7 G [0,1 ). For a smooth step, 5 is defined 
by the ratio h/d = h/d e [0,oo) and the formula (2.13) can be explained as follows 


7(X) = 


r 1 , h/d- 

\ 0, h/d- 


oc for h 7 ^ 0 
0 


, for f{X) = fs{X,h),X G [-d/2, d/2]. (2.14) 


When 7=1 the step is a sharp whereas for 7 = 0 the smooth step tends to a flat plate. 
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Figure 2. Streamwise velocity profiles at 3 different streamwise locations: (a) Res* = 821, (b) 
Res* = 866 and (c) Res* = 897. The physical parameters corresponding to each case are from 
Case E in Table 


3. Numerical approach and results 

3.1. Numerical strategy 

A spectral//ip element discretisation, implemented in the Nektar++ package, is used in 
this work to solve the linear as well as nonlinear Navier-Stokes equations. A stiffly stable 
splitting scheme is adopted which decouples the velocity and pressure fields and time 


integration is achieved by a second-order accurate implicit-explicit scheme (Karniadakis 
et Cantwell et a/. [[2^5 ). 

For 2D simulations, a convergence study by p-type refinement is performed to demon¬ 
strate mesh independence is achieved throughout this study. For 3D calculations, pro¬ 
vided the same 2D mesh in the x — y plane is used with the addition of a hybrid Fourier- 
Spectral/Zip discretisation in the third direction yielding a hybrid Fourier/spectral/Zip 
discretisation of the full 3D incompressible Navier-Stokes equations. 

To obtain the best fidelity in representing the curved surface of the smooth step, high- 
order curved elements are defined by means of an analytical mapping. The governing 
equations are then discretised in each curved element by seventh-order polynomials. The 
choice of polynomial order is guaranteed by the mesh independence study. Figuresand 
compare the horizontal and vertical velocity profiles over a smooth step in Res* = 866 are 
given at different positions for polynomial order ranging from six to eight. In the whole 
domain, the relative error of velocity fields is lower than 10“^, which is consistent 
with the convergence tolerance of the base flow generation defined by 

\\{d^u,d^v)\\J\\{u,v)\\o ■ < 10“®, 

where || • ||o means the standard norm, 9/ denotes the discrete temporal derivative 
and Tc is the convective time scale. Once 2D steady base flows are generated by using the 
non-linear Navier-Stokes equations (NSEs), the TS waves are simulated by the linearised 
Navier-Stokes equations (LNSEs). As discussed below, for base flow generation, the inlet 
position is located sufficiently far from the first step to allow the base flows to recover 
the Blasius profile. Following the experimental methodology used by [Downs fc Fransson 
(2014), the TS waves are excited by periodic suction and blowing on wall. 
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v/ 


'^max,modes=8 


Figure 3. Wall-normal velocity profiles at 3 different streamwise locations: (a) Re§* — 821, (b) 
Res* — 866 and (c) Res* — 897. The physical parameters corresponding to each case are from 
Case E in Table 


3.2. Linear analysis in a narrow unstable regime of the neutral stability eurve 

In order to detect the effect of smooth steps on the stability of the boundary layer, 
three different high frequency perturbations E (g{ 140, 150 160}) are excited within 
the boundary, up-stream of the unstable region, by periodic blowing and suction (see 
figure]^. Hereafter, for each perturbation frequency studied, let Aq indicate a reference 
maximum TS mode amplitude at the lower branch of the neutral stability curve in a 
fiat plate boundary layer. The contours of Iftl/Ho are given in figure]^ for three different 
non-dimensional frequencies {E G {140,150,160}) for four different single smooth steps of 
different height and smoothness all located in Re^^i = 680. A summary of the parameters 
of these computations can be found in table 

First, we observe that around a smooth step, the TS mode is energised and subse¬ 
quently weakened (figure IHal -d3)). By energising we mean that around the step, there 
exists a local maximum of |fi|/Ao. Secondly, a higher step height h gives a stronger local 
maximum. Thirdly, increasing the height of the step moves the location of maximum 
amplitude of the TS mode downstream. For example for the an excitation frequency of 
E = 150 the maximum is located at Res* = 920 for h = 5.48% (figure (^ 2 )) and at 
Res* = 950 for h = 30% (figure]^ (d 2 )). 

The result of our linear analysis reveal distinct behavior from the results for a sharp 


step (7 = 1) of height = 0.235 {h = 8.24%) given by 


Worner et al. 


(2003) where 


they claim that a forward facing sharp step showed a stabilising effect without a local 
destablisation regime despite a separation bubble being reported in front of the step. 
Recently, through further numerical calculations, Edelmann & Rist (2015) observed that, 
for subsonic Mach numbers and large height scale {h/S* = 0.94 or h = 33%), two 
separation bubbles, in front and on top of the step, were observed and strong amplification 
of the disturbances was found in front of and behind the step. The results from the linear 
analysis and, more generally the results from the DNS presented in §3.4, underline the 
attractiveness of replacing a forward facing sharp step (7 = 1 with large h hy a. smooth 
one 7 < 1 because the smooth step with the same h does not lead to a separation bubble. 
We attribute the discrepancy between the current findings and the work of [Edelmann fc: 


Rist (2015) to the step-induced separation bubble that has strong destabilising effect on 


the TS mode. Furthermore, for a forward-facing smooth step of fixed height, the position 
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of global maximum value of a coutour varies with respect to the giveu frequeucy, which is 
cousisteut with the relative positiou of the step with respect to the positiou of the upper 
brauch wheu chaugiug frequeucy iu the ueutral curve diagram ([^cl-c3) for example). 
For two isolated smooth steps, a similar pheuomeuou is observed iu figure [^cl-c3) except 
that two distiuct local maxima are observed arouud each smooth step for large h. From 
figures aud it cau be coucluded that both for siugle- aud two-step coufiguratiou 
global maximum values of \u\/Ao depeuded ou frequeucies, h aud smoothuess. 

By extractiug maximum values of the euvelope of the TS mode at each streamwise 
locatiou, the Aj profiles for each case are showu iu figure We uotice that wheu 
height scale of oue siugle smooth step is uot high euough, it does siguificautly euergise 
the TS wave aud the euergised regiou is limited to the viciuity of the the step. From 
figure [^ai -as), we observe that wheu h < 20 %, the TS waves are stabilised dowustream 
for small h. Iu figure a comparisou of the TS modes at three differeut locatious are 
showu, further elucidatiug the stabilisatiou effect. We hud that for a siugle smooth step, 
wheu frequeucies are chauged, behaviours of the TS modes are differeut from each other. 
Wheu h < 20%, for all three frequeucies, the global maximum value of the TS mode 
does uot exceed the global maximum values of the TS modes iu the correspoudiug fiat 
plate bouudary layers. Aualysiug the frequeucies iu figure wheu h ^ 20%, the most 
amplified TS modes correspoud to the lowest frequeucy T = 140. The reasou for this 
lies iu that for a fixed step positiou, choosiug a low frequeucy iudicates that the step is 
moved towards the lower brauch of the ueutral stability curve of a flat plate bouudary 
layer. This is ouly a heuristic explauatiou, siuce large step heights have a siguificaut 
impact ou the positiou of the lower brauch, which cau lead to deformatiou of the shape 
of the ueutral stability curve aud moviug the ueutral positiou upstream with respect to 
that of a flat plate bouudary. This uou-local pheuomeuou cau be seeu from a comparisou 
betweeu the TS mode iu a flat plate bouudary aud the TS mode iu a bouudary layer 
over a large height scale smooth step as illustrated iu figure [^ai -as). 

From figure [^ai- as), assumiug the local growth (or destabilisatiou) of the TS wave 
does uot trigger auy uouliuear pheuomeua, a smooth step with a low height scale is 
uot harmful for the TS wave with a high frequeucy (140 ^ ^ 160). Iu fact, to some 

exteut, a bouudary layer cau beuefit from a smooth step siuce the uet iustability cau be 
reduced. Iu figure [^bi-bs), with the same frequeucies, euvelopes of the TS waves over 
two isolated smooth steps are showu. The positiou of the two smooth steps cau be fouud 
from Table as schematically illustrated iu figure |^b). For the frequeucies cousidered 
here, the secoud steps still lie iu the uustable regime of the ueutral stability curve. We 
observe that, surprisiugly, the secoud steps do uot locally lead to further amplificatiou 
of the TS waves wheu h < 20 %; iu coutrast, the amplitudes of the TS waves are further 
damped compared to the amplitude of the TS waves amplitudes over both the siugle 
step bouudary layers aud a fiat plate bouudary. A destabilisiug of the TS mode is ouly 
iutroduced by large height scale steps. 

From the results of the liuear stability aualysis preseuted above, we deduce that stabili¬ 
satiou or destabilisatiou behaviours of smooth steps strougly depeuds ou the smoothuess 
parameter 7 as well as the height. Iu a suitable rauge of parameters the iufiueuce of 
smooth steps ou a bouudary layer cau be domiuated by stabilisatiou effect. To couclude, 
iu the uustable regime of the ueutral stability curve, which is typically uear the leadiug 
edge, if a step-like structure is iuevitable, a smooth step will have less of destabilisiug 
effect ou the TS waves thau a sharp oue. 
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Res* Res* 


Figure 4. Positions of exciters (+) and one single smooth steps (a) and two smooth isolated 
steps (b). Three different frequencies corresponding to steps’ positions are given by horizontal 
dashed lines. 


Case 

Res* 

Res,^ 


El 

E2 

Es 

h% d 

7 X 10^ 

LxISqq 

Ly/699 

A 

320 

680 

786 

140 

150 

160 

5.48 4 

0.74 

250 

30 

B 

— 

— 

— 

— 

— 

— 

10.96 — 

2.99 

— 

— 

C 

— 

— 

— 

— 

— 

— 

20.00 — 

9.97 

— 

— 

D 

— 

— 

— 

— 

— 

— 

30.00 — 

22.44 

— 

— 

Table 1. Parameters for smooth steps where Res^. 

, Res* and Re5* 

Cl C2 

are, respectively, the inlet 


Reynolds number, the Reynolds number at the centre of the first step and the Reynolds number 
at the centre of the second step . T denotes the non-dimensional perturbation frequency. Lx and 
Ly denote streamwise extent and height of the domain for which the 2D base flow field obtained 
was independent of domain size. 


3.3. Linear analysis of smooth steps at lower exeitation frequeneies 

In the previous section, the investigations of the effect of a smooth forward facing step 
on the stability of a boundary layer focused on the high frequencies T G {140,150,160}. 
The unstable regions corresponding to these frequencies are relatively narrow compared 
with lower frequencies. For example for T = 160 the unstable region of the TS-mode 
on the flat plate ranges from Res^ = 580 to Res* = 830, whereas for a perturbation 
frequency oi T = 100 the unstable region ranges from Res* = 700 to Res* = 1100. In 
this section, we consider the TS wave with a frequency = 100 in a boundary layer over 
a single smooth step and the cases with the parameters given in Table Physically, the 
size of steps at is kept the same as that of steps at Res^ • We are now interested in 

assessing the effect of the position of the step on the excitation of the TS mode. In figure 
9, the comparison of the TS envelopes in the boundary layers is given. Clearly, figure 
9 ^a) indicates that small h cannot induce significant destabilisation. If h is increased as 
seen in figure [^b)-(d), a destabilisation effect emerges and larger h give rise to larger 
global maximum amplitude of the TS wave. Meanwhile, figures [^c) -(d) clearly indicate 
that moving the location of the smooth step downstream, far from the upper branch of 
the neutral stability curve, it induces more amplification than that of a more upstream 
location. 

This means that when the TS wave has a low frequency, smooth step could play a 
destabilisation role for the TS wave. From figure we learn that growth rate of the 
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Figure 5. Contour plots of Iftl/Ao and the physical mrameters corresponding to (a^), (b^), 
(c^) and (d^) are from Case A, B, C and D in Tabled where # denotes the number 1, 2 or 3, 
which corresponds to frequency The step position is determined by Re^c^. The grey cross 
+ indicates the location of the maximum amplitude of the TS-wave and the vertical grey line 
represents the location of the forward facing smooth step. 
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Figure 6. Contour plots of Iftl/Ao and the physical parameters corresponding to (a^), (b^), 
(c^) and (d^) are from Case A, B, C and D in Tabled where # denotes the number 1, 2 or 3, 
which corresponds to frequency The step positions are determined by Re^<^\ and Re^<^^ . The 
grey cross + indicates the location of the maximum amplitude of the TS-wave and the vertical 
grey line represents the location of the forward facing smooth step. 
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Figure 7. The physical parameters corresponding to A, B,C and D are from Case A, B, C 
and D in Table [b (a^) and (b^), respectively, denote single step and two steps where where ^ 
denotes the number 1, 2 or 3, which corresponds to frequency in Table The vertical grey 
lines represents the location of the forward facing smooth step. 


TS waves in a flat boundary at Res* is greater than that at Res*. We attribute the 
greater destabilisation effect induced by the step at Res* to the greater growth rate at 
that position in a flat plate boundary layer. This phenomenon is further considered in 
figure 


10 


for = 20 and 30 for the same nondimensional scales h and d at Res* = 866 
and Res*^ = 988. The correlation between the position of the step and the TS mode 


amplification acts as a guideline for choosing the location of a smooth step. For large /i, 
the effect of a step located in a larger growth rate region gives rise to larger amplification 
of the TS wave. However, it is worthy to note that when height scale is reduced, there 
no longer exists a strong destabilisation effect from a single smooth step, regardless of 
its position within the unstable region. 

To compare the effect of one and two steps additional simulations are carried out for 
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Figure 8. Comparison of the TS modes at different location over a single smooth step. The 
parameters corresponding to A, B, C and D in the figures are from Case A, B, C and D in 
Tableand T — 150. (Left-Top) and (Left-Bottom): Re§ = 680; (Mid-Top) and (Mid-Bottom): 
Re§ = 750; (Right-Top) and (Right-Bottom): Re§ = 800. u and v are, receptively, normalised 
by max('ifc/) and max('5/) from the corresponding flat plat boundary. 


two smooth steps located in regions of large growth rate of the TS wave. The amplification 
A/Aq is plotted as as function of Res* in figure [T^a) for the single step an d[^b) for the 
two steps. We observe that for small h{< 20%), the second steps do not locally induce 
more destabilisation effect than the single step. Case A in figure [T^b) indicates that for 
small h{< 20%), there exists a week stabilisation effect over the second step. However, a 
strong destabilisation effect is observed from h{'^ 20%) corresponding to cases C and D 
in figure [T^a-b) . Both for the single and the two step case, increasing the height of the 
step moves the location of the maximum amplitude of the TS mode to higher Res *. 


3.4. Does a smooth step signifieantly amplify the low-frequeney TS waves ? 

The studies by Worner et a/.| (2003); Edelmann & Rist (2015) considered the effect of low 
frequency TS modes for a long range from the leading edge of the flat plate with Res^ 
(i.e. Res^ = 2200 1500). In order to further assess the effect of a forward facing smooth 

step on the amplification of the TS wave at high Res^ regime where the N factor is close 
to 8, we consider a TS wave with = 42. The excited TS wave is further amplified by 
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Case 

Res* 

1 


R 

h% 

d 

7 X 10^ Lr/d99 

Ly/Sgg 

A 

388 

866 

100 

5.16 

4 

0.66 

287 

30 

B 

— 

— 

— 

10.32 

— 

2.65 — 

— 

C 

— 

— 

— 

20.00 

— 

9.97 — 

— 

D 

— 

— 

— 

30.00 

— 

22.44 — 

— 


Res* 

Res*, 

R 

h'(=M§9/4)% d' 

= dSgg/ (^§9 

7 X 10^ Lr/d99 

Ly/Sgg 

A' 

388 

988 

100 

4.52 

3.51 

0.66 

287 

30 

B' 

— 

— 

— 

9.05 

— 

2.65 — 

— 

C' 

— 

— 

— 

17.54 

— 

9.97 — 

— 

D' 

— 

— 

— 

26.31 

— 

22.44 — 

— 


Table 2. Parameters for smooth steps: Res* is the inlet Reynolds number, Res* and Res*^ 
indicate two different locations with respect to two same-size single steps. T denotes the non-di¬ 
mensional perturbation frequency. 



Res* Res* 



Figure 9. Comparison for steps at different locations and the parameters for (a)-(d) are given 
in Table 1^ (^-D) one single step at Res*] (A'-D') one single step at Res*, - The vertical grey 
line represents the location of the forward facing smooth step. 


the recirculation bubble inside the indentation, located in Res* = 1519, far upstream of 
the transition criterion position of N = S. We then add a smooth forward facing step 
downstream of the indentation in the unsteady region of the TS mode. This configuration 
has its practical interest in roughness prompted transition. The indentation is defined by 


fr 


-hr • COS (nX^/xf , X’’ e [-A/2, A/2], 

0, X-^[-X/2,X/2], 


(3.1) 





















16 


H. Xu, J.-E. W. Lombard, S. J. Sherwin 


Case 

Res* 

l 


Res-^ 

E 

h% 

d 

7 X 10^ 

Lx / bgg 

Ly/Sgg 

A 

388 

988 

1096 

100 

5.00 

4 

0.58 

287 

30 

B 

— 

— 

— 

— 

10.00 

— 

2.30 

— 

— 

C 

— 

— 

— 

— 

20.00 

— 

9.26 

— 

— 

D 

— 

— 

— 

— 

30.00 

— 

20.84 

— 

— 


Table 3. Parameters for smooth steps where Res*, R^s*^ and Res*^ are, respectively, the inlet 
Reynolds number, the Reynolds number at the centre of the first step and the Reynolds number 
at the centre of the second step. T denotes the non-dimensional perturbation frequency. 




Figure 10. The physical parameters corresponding to A, B,C and D are from Case A, B, C 
and D in Table (a) one single step at ; (b) two steps at Re^<^\ and R^sl^- vertical 

grey line represents the location of the forward facing smooth step. 


where X'^ = x — The parameters of this computation are given in Table 

Note that the width scale A is comparable with the corresponding TS wavelength. A 
separation bubble is induced in the indentation region and when a base flow undergoes 
this distortion, the TS wave is strongly amplified. In figure miL eft), we observe that 
the TS waves are strongly amplified around the indentation and the N-factor in the 
computational domain reaches N= 8 . Downstream of the roughness, for Res^ > 1800, the 
smooth forward facing step only has weak, local, stabilising or destabilising effect for all 
three cases (see figure [TT]^ Right)). 

For the case /i = 5, the envelope of the TS wave is not significantly modified by the 
smooth step. Without talking about any possible stabilisation effect, destabilisation effect 
for such a low frequency is negligible. That is to say, for a TS wave with a low frequency 

= 42, a smooth step does not significantly impact the N-factor. In contrast, a strong 
destabilisation of the TS wave by the sharp step (7 = 1 ) for low frequency is reported 
by Edelmann & Rist (2015) where the separation bubble induces a strong increase in 
N-factor from N = 4 to N = 6 because of the existence of separation bubbles in front 
of the step. This contrasts with the weak destabilisation influence of a smooth step on 
the TS wave in a boundary layer where negligible variation of the N-factor can be seen 


in figure 11 even for a large h > 20 % step height. 


3.5. DNS investigation of the effeet on K- and H-type transition of two forward faeing 
smooth steps exeited by a high frequeney TS wave 

The TS mode is two dimensional, which is why linear stability analysis of the boundary 
layer can be conducted in two dimension. However secondary instabilities and transition 
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Case 

Res* Res* 

Res* 

R 

hr% 

h% 

A 

d 

hx !^99 

Ly /899 

o 

596 1519 

1885 

42 

74.77 

0 

5.5 

4 

312 

30 

A 

— — 

— 

— 

— 

5.00 

— 

— 

— 

— 

B 

— — 

— 

— 

— 

10.00 

— 

— 

— 

— 

C 

— — 

— 

— 

— 

15.00 

— 

— 

— 

— 

D 

— — 

— 

— 

— 

30.00 

— 

— 

— 

— 


Table 4. Parameters for a wall with an indentation and a smooth step. and Res* are, 

respectively, the inlet Reynolds number, the Reynolds number at the centre of the indentation 
and the Reynolds number at the centre of the smooth step. T de note s the non-dimensional 
perturbation frequency, hr and A are used to define the indentation (3.1), which are normalised 
by the boundary layer thickness at the centre position of the roughness. 




Res* 


Res* 


Figure 11. Effect of a smooth step for the TS wave with a low frequency. N — ln(T/Ao). 
The roughness for amplifying the TS wave is located at the position of the first vertical line 
and the step is located at the position of the second vertical line. (Left) Overview of the TS 
waves’ envelopes; (Right) Local view of the TS waves’ envelopes around the smooth steps. The 
parameters are from Table The vertical grey line represents the location of the forward facing 
smooth step. 


to fully developed flow is a highly three dimensional. To this end further investigation 
of the influence of smooth steps on two transition scenarios is achieved with a hybrid 
Fourier-Spectral//ip discretisation is employed to solve 3D Incompressible Navier-Stokes 
equations. The spanwise direction was assumed to be periodic and discretised by 80 
Fourier modes and the streamwise and wall normal plane was discretized using 5576 
elements (quad and triangle) within which a polynomial expansion of degree 7 is imposed. 
K- and H-type transitions are simulated for the flow over a flat plate and with smooth 
steps. A Blasius profile is imposed at the inflow and, for both scenarios, a wall-normal 
velocity along the disturbance strip is prescribed by the blowing and suction boundary 
condition ( Huai e^'^|1997 ) , 


v{x, zfl) = A - f{x) • sm{ujAt) + B • f{x) • g{z) • sm{ujBt), 


(3.2) 


where uja and ujb are the frequencies of the 2D TS wave and the oblique waves, re¬ 
spectively. A and B are the disturbance amplitudes of the fundamental and the oblique 
waves. The function f{x) is defined by ( Fasel fc Konzelmann||199Q' ) 


f{x) = 15 . 1875 ^^ - 35 . 4375 ^"^ + 20 . 25 ^^ 


( 3 . 3 ) 
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Res* 

% 

Resi^ 


Ea 

Eb 


B/C/co 

h% 

Lx / dgg 

Ly / (igg 

Lzjdgg 

y^zjdgg 

T 

320 

680 

786 

150 

150 

75 

150 

0.5% 

0.03% 

0 

5.48 

340 

30 

8 

4 

20 

— 

— 

— 

— 

75 

150 

75 

— 

— 

5.48 

12.79 

12.79 

— 

— 

— 

— 

— 


Table 5. Parameters used for the DNS simulations. Ea and Eb denote non-dimensional per¬ 
turbation frequencies of the disturbance strip. Ajlloo and Bjlloo are the relative amplitudes of 
the disturbance amplitudes of the fundamental and oblique waves, respectively. The spanwise 
Lz extent of the domain is expressed as function of the boundary layer thickness Sgg. T is the 
finial nondimensional time length scale which we simulate, which is defined by T = tUoo/L. 



Figure 12. Overview of the computational setup with the Blasius boundary layer profile at 
the inflow, the disturbance strip and two smooth steps used for the DNS. 


with the parameter ^ 




" X — Xi 
X2 - X 

^ X2 X^ui 


for Xi ^ X ^ Xm 
for Xm ^ X ^ X2 


(3.4) 


where Xm = {xi + X2)/2 and g{z) = co^{ 2 t:z / \z) with the spanwise wavelength A;^. 
At xi and X 2 , and Res* are equal to 591.37 and 608.51, respectively. For K- 

type transition, the oblique waves have the same frequency as the two dimensional wave 
{uja = 00b)- Also, for H-type transition, the oblique waves are sub-harmonic {ojb = ojaI‘2)- 
All parameters used in the investigation are given in Table Note that the frequency 
used for the perturbation are consistent with the parameters used for linear analysis 
in Table The schematic of the computational domain is illustrated in figure It is 
worth mentioning that, in order to guarantee that nonlinear calculations are converged, 
all unsteady simulations terminate when the nondimensional time scale T is equal to 20. 


3.5.1. Effect of smooth forward facing steps on K- and H-type transitions 

The validation of the DNS results for both K- and H-type transitions is corroborated 
by recovering the aligned arrangement of the A vortices for the K-type transition (see 
figure(a)) and staggered arrangement for the H-type transition (see figure (b)) as 
experimentally observed by Berlin et al. (1999) for the fiat plate boundary layer. Figure p^ 
shows the evolution of the skin-friction coefficient versus Res* for two different normalized 
height scales. We observe that, for a fiat plate boundary layer, the skin friction coefficient 
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Figure 13. Instantaneous contours of stream-wise velocity in xz-p\?me at height y — 0.6(599 in 
Res* G [963,1111] for the K- (a) and H-type (b) transition scenarios for a flat plate h = 0. 


diverges from that of the Blasius boundary layer at the location where A vortices appear, 
as illustrated in figure 13 The streamwise evolution of the skin-friction (see figure 14) 
shows the K-type transition is fully inhibited by the two smooth steps whereas the H- 
type transition is delayed. Additionally, increasing the height h (< 20 %) further reduces 
the skin friction coefficient Cf in both scenarios. The observation of these phenomena 
supports the result of linear analysis. 

To gain further insight into the different impact on two transition scenarios of the two 
forward facing smooth steps we consider the energy growth of the main modes. We label 
these modes using the notation (cj,/3) ( [Berlin et al. 1999), where uj and (3 are respec¬ 
tively, the frequency and spanwise wavenumber each normalized by the corresponding 
fundamental frequency/wavenumber. It has been observed that the K-type transition 
scenario has the main initial energy in the (1,0) mode. The (1, ±1) mode also generates 


the (0, ±2) mode with a small amplitude through non-linear interaction (Berlin et al. 
1999| ). At the late stage, the ( 0 , ± 2 ) mode can grow to an amplitude comparable to that 


of the (0,±1) mode. Laurien & Kleiser (1989) and Berlin et al. (1999) have shown that 
the initial conditions for the H-type transition has the main energy in the (1,0) mode 
with a small amount in the oblique subharmonic (1/2, ±1) mode. The important mode 
is the vortex-streak ( 0 , ± 2 ) mode, which is nonlinearly generated by the subharmonic 
mode and vital in the transition process. As illustrated in figure [T^ the (0, 2) mode plays 
a significant role in the late stages of transition for both transition scenarios with the 
two smooth steps. For the K-type scenario, in the transition regime, the energy of mode 
(0, 2) grows and exceeds the energy of mode (0,1). For the fiat plate, the energy of mode 
(0,1) finally grows again until turbulence occurs. The energy of mode (0, 2) with the two 
smooth steps is less than that of mode ( 0 , 2 ) for the fiat-plate and increasing normalised 
smooth step height h yields stronger reduction of the energy. A similar reduction in en¬ 
ergy is also observed for mode (0,1). Furthermore, the energy of mode ( 0 , 2 ) exceeds that 
of mode ( 0 , 1 ) in Res* = 1080 for h = 5.48% to Res* = 1100 for h = 12.79% (figure 15 


(a)) and from this points onwards the energy of mode (0,1) decays. Based on the results 
presented in figure (a), we observe that the spanwise modulation induced by mode 
( 0 , 1 ) with energy decaying on the smooth steps leads to the stabilisation of the boundary 
layer. 


4. Conclusion 

In this paper we introduce the smoothness of a step 7 and show that it affects on 
destabilising effect of step and on an otherwise fiat plate boundary layer. Linear stability 
analyses are conducted with high 140 < < 160, medium = 100 and low iF = 42 fre¬ 

quency forcing with respect to the neutral stability curve of the flat plate boundary layer. 
One and two step configurations with different heights and smoothness are considered. 
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Figure 14. Comparison of time- and spanwise-averaged skin friction versus streamwise position 
Rex for K- (a) and H-type (b) transition scenarios for a flat plate (—) and two of height 

h = 5.48% (-) and h = 12.79% (-). The skin-friction profile of the Blasius boundary layer 

(•) is given for reference. 




Res* 


Res* 


Figure 15. Comparison of the energy in modes (0,1) and (0,2) versus streamwise position Res* 

over a flat plate (—) and two smooth steps with height h = 5.48% (-) and h = 12.79% (-) 

for K- (a) and H-type (b) transition scenario. 


For the low frequency forcing case an indentation is added upstream of the forward fac¬ 
ing smooth step and the effect on the stability of this highly perturbed flow is reported. 
Finally direct numerical simulation of K- and H-type transition scenarios are undertaken 
for a high-frequency forcing case to confirm the results from the linear analysis. 

The net effect of a smooth forward facing step on the stability of the TS-mode depends 
on height, width and flatness. From an application point of view, small height smooth 
steps are safe for TS waves of frequency ranging from T G [42,160]. For h < 20% both 
the single and two forward facing smooth step configurations lead to a stabilising effect of 
the TS-mode of up to 5% and 10% respectively. For h > 20% the smooth step 7 < 1 has 
a weaker destabilising effect than previous reports show for a sharp step 7 = 1 . Further 
investigation indicates that for the TS wave with a low frequency, small height smooth 
steps can be safe and don’t significantly amplify the TS wave. This result contrasts with 
that of a sharp forward facing step and we attribute the large amplification of the TS 
wave by a sharp step to separation bubbles. This effect is similar to the destabilisation 
effect induced by the separation bubbles in an indentation or behind a hump or a bump 

From a global stability point of 


(Gao et al. 2011 Park & Park 2013 


Xu et al. 


2016). 
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view, a smooth step, even with a large height, cannot lead to a destabilisation in front 
of or over the step because of the absence of recirculation bubbles meaning no global 
instability can be introduced. The results obtained by DNS, for a high-frequency forcing 
T = 150, support the conclusion that smooth steps can have non-negligible and positive 
impact on the stability of the boundary layer. For h = 5.48 and h = 12.79 the transition 
to fully developed turbulent state is even delayed for the H-type transition scenario and 
suppressed for the K-type scenario. 

More generally these results suggest a smooth forward facing step can have a significant 
influence on the stability of the boundary layer. 
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